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Abstract 

Using techniques from hopping expansion we identically map the 
lattice Schwinger model with Wilson fermions to a model of oriented 
loops on the lattice. This is done by first computing the explicit form 
of the fermion determinant in the external field. Subsequent integra- 
tion of the gauge fields renders a sum over all loop configurations with 
simple Gaussian weights depending on the number of plaquettes en- 
closed by the loops. In our new representation vacuum expectation 
values of local fermionic operators (scalars, vectors) can be computed 
by simply counting the loop flow through the sites (links) supporting 
the scalars (vectors). The strong coupling limit, possible applications 
of our methods to 4-D models and the introduction of a chemical po- 
tential are discussed. 
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1 Introduction 



The evaluation of the fermion determinant in lattice field theories poses 
a highly non-trivial problem. In particular for an odd number of flavors, 
or when a chemical potential is added, no fast numerical algorithms are 
known. As a possible way out several studies of alternative representations 
of the fermion determinant in terms of fermion loops can be found in the 
literature |l]]-[10]. Loop representations are not only useful as a tool for 
numerical simulations, but also provide interesting insights into the nature of 
the fermion determinant, in particular they provide a solution of the fermion 
sign problem (compare []Tl|,|l^1 for recent developments). At least for abelian 
gauge fields it is furthermore possible to perform a duality transformation 
for the link variables and in 2 dimensions they can even be integrated in 
closed form. The resulting representation in 2-D is a discrete sum over filled 
Wilson loops. 

The historic approach to finding a loop representation for the fermion 
determinant is to first expand the Boltzmann factor for the fermions. Since 
the fermions are represented by Grassmann variables the power series for 
the Boltzmann factor has only finitely many terms. In a second step the 
Grassmann integral at each site is saturated with all possible combinations 
of the terms from the series for the Boltzmann factor. The result can be 
viewed as a sum over closed loops of links. We remark, that for gauge 
theories products of link variables arranged along closed loops are the only 
possible gauge invariant terms and thus structures of closed loops can be 
expected already from gauge invar iance. 

The local saturation of the Grassmann integral is however a tedious 
algebraic problem. Thus most of the results were obtained for staggered 
fermions where the additional structure due to spinorial degrees of freedom 
is avoided fl]]-f|. The resulting model however still describes more than 
one fermion due to the doubling problem, which is only reduced but not 
completely eliminated by the staggered formulation. 

The first result for Wilson fermions is Salmhofer's mapping of the strongly 
coupled Schwinger model to a 7 vertex model [0]. The techniques used there 
were later generalized by Scharnhorst and applied to the strongly coupled 
Schwinger model with 2 flavors and the 2-D lattice Thirring model [|]. 
Thus for purely fermionic models in 2-D also loop representations for Wilson 
fermions have been found for some models. However, no loop representation 
for the fermion determinant in an external field has been found for Wilson 
fermions while for the staggered case such results are known 0. 
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In a recent publication we have attacked this problem using a new ap- 
proach [0, [nj. The central idea is to first perform a hopping expansion of 
the fermion determinant (see e.g. |13| for an introduction and [14|] for an ap- 
plication to the Schwinger model). The hopping expansion essentially gives 
a loop representation for the free energy of fermions in an external field. 
Thus it already provides the natural variables for the loop representation. 
The remaining problem is to expand the Boltzmann factor in order to bring 



down the loops from the exponent. In |], 1C] this method has been applied 
to the case of 2-D Wilson fermions interacting with an external scalar field 
through a Yukawa coupling. 



In this article we treat the case of the lattice Schwinger model with 
Wilson fermions using the hopping expansion technique. The fermion de- 
terminant in the external U(l) gauge field will be written as a sum over 
oriented loops dressed with the gauge variables sitting on the links of the 
loops. In this representation we then integrate explicitly the gauge fields 
and the partition function turns into a discrete sum over filled Wilson loops. 
Vacuum expectation values of Wilson loop type observables can easily be 
taken into account by adding the Wilson loops to the loops from the deter- 
minant before integrating the gauge fields. VEV's of local fermionic observ- 
ables (scalars and vectors) also turn out to have rather simple expressions 
in the loop formalism. They can be obtained by differentiating the fermion 
determinant (before performing the gauge field integration) with respect to 
the mass term for the scalars or with respect to the link variables for the 
vectors. In the loop picture these observables then correspond to the loop 
flow through sites (for the scalars) and links (vectors). Thus in a numerical 
simulation of the loop model the evaluation of n-point functions for these 
observables becomes exceedingly simple. Finally we remark, that a chemical 
potential can be included in a natural way by giving the pieces of loops run- 
ning forward in time a weight different from the weight for pieces running 
backward in time. 

The article is organized as follows: In the next section we rederive the 
hopping expansion of the fermion determinant in a form which is convenient 
for our later manipulations. In Section 3 we expand the Boltzmann factor 
and construct the loop representation of the determinant. This is done by 
first analyzing several examples for different contributions and subsequently 
the general rules of the loop calculus are discussed. In the end of Section 
3 we also briefly address the role of a chemical potential. In Section 4 we 
discuss the relation of our new representation to known results which can 
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be obtained as limiting cases of our formulas. In particular the coupling 
to an external Yukawa field and the strongly coupled Schwinger model are 
considered. In Section 5 we integrate out the gauge fields and give the final 
representation for the partition function of the lattice Schwinger model with 
Wilson fermions in terms of filled loops. In Section 6 we discuss the loop 
representation of some observables and the article ends with a discussion in 
Section 7. 



2 Setting and hopping expansion 

We study the lattice Schwinger model with Wilson fermions. In addition to 
the gauge fields we also couple a real scalar field 9 to the fermions through 
a Yukawa-type interaction. The scalar field 9 will serve as a tracking device 
for the algebraic manipulations below and also allows for a consistency check 
by comparing our new representation with the known loop representation 



for the Yukawa- type coupling |10|. In addition we will use 9 as a source 
term for generating n-point functions. 

The Wilson lattice action is a bilinear form S = ^2 x ^ y il)(x)K(x : y)4>{y) 
with kernel 



K(x,y) = 2 + m + 9(x) 



->x,y 



±2 

E 

m=±i 



J x+jl,y i 



where we defined 



r ±M = -[it^I 



1,2 . 



Here a±, o"2 are Pauli matrices. U^(x) £ U(l) are the link variables supported 
on the links (x, [i] of the lattice, with U-^(x) = U^(x — jl)* , (fj, = 1, 2), i.e. for 
hopping in backward direction the link variables are complex conjugate. The 
sum in the bilinear form runs over the whole lattice, x, y G A. For simplicity 
we assume that our lattice A is a finite rectangular piece of 2Z 2 with open 
boundary conditions, i.e. hopping terms that would lead to the outside of 
our lattice are omitted. We remark that all our results also hold for A 
being a cylinder with a compact time direction with anti-periodic boundary 
conditions and open boundary conditions for the space direction. Such a 
geometry is interesting for analyzing the model at finite temperature. Let's 
define 

M(x) = 2 + m + 9(x) , 
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and assume that 9(x) is such that M(x) 7^ for all lattice points x. This 
is a purely technical assumption due to the particular techniques we use for 
computing the determinant. The final result will be a finite polynomial in 
the 9{x) and the above restriction is irrelevant then. We now can write for 
the fermion determinant 



/ dipd^exp — y~] ip(x)K(x,y)ip(y) = detK 

( OO -i \ 
-J2-^H n \ . (2) 
n=l n / 

In the last step the hopping expansion was performed, i.e. the determinant 
was expressed using the well known trace-logarithm formula and the loga- 
rithm was expanded in a power series. The hopping matrix H is defined 
as 

H(x,y) = E r p||^r ( 3 ) 
The series in the exponent of @ converges for \\H\\ < 1, and the whole 



expression represents the determinant then (see e.g. [lq] ). In 16] it is shown 
for the case of 9(x) = 0, that \\H\\ < 1 for m > 0. It is straightforward to 
generalize this argument and prove that for our case we have < 1 for 
m + 9{x) > 0. To summarize: the series in (2) converges for m + 9{x) > 
and the right hand side of (|2|) represents the determinant for this range of 
parameters. 

Due to the Kronecker delta in @, the contributions to TrH n are sup- 
ported on closed loops on the lattice, and since closed loops are of even 
length (for the cylinder geometry discussed above we assume an even length 
for the compact direction), the contributions for odd n vanish. For even 
n = 2k we obtain 

T^ 2fc = e e n ^iir.. (4) 

Here is the set of all closed, connected loops of length 2k and base 

point x. In the first product we collect all the link variables U fJi (y)/M(y) on 
the links visited by the loop according to their orientation. The last term in 
(||) is the trace of the ordered product of the hopping generators T v as they 
appear along the loop I. We remark, that T± U T TU = 0, which implies that 



4 



whenever a loop turns around at a site and runs back along its last link this 
contribution vanishes. Thus all these back-tracking loops can be excluded 
from C^ k \ 

Evaluating the trace over the matrices T u for a given loop is the re- 
maining problem in the hopping expansion. It first has been solved for the 
infinite lattice in [lTj by realizing that the Pauli matrices give rise to a rep- 
resentation of discrete rotations on the lattice. Alternatively one can build 
up the loops using four basic steps and compute the trace in an inductive 



procedure along the lines of [18]. The latter method can be extended to the 
case of lattices with compact dimensions. When choosing a cylinder with 
anti-periodic boundary conditions in time direction and open b.c. for the 
space direction the result equals the result for the rectangular lattice with 
open boundary conditions. Thus the result for both of the discussed sets of 
boundary conditions is 

Trll^ = -(-ir«(^) CW . (5) 

By s(l) we denote the number of self-intersections of the loop I and c(l) 
gives its number of corners. We remark, that for e.g. a torus additional sign 
factors for topologically non-trivial loops would appear. The result (0) is 
independent of the orientation of the loop. Inserting (|j) and (]5|) in (|2|) we 
obtain 



det K = Yl M(x) 2 exp 



xeA 



EiEE (-^f n W> 

k=i 2k v eA ie£ (2*) V V 2' Mel M (z) 

(-iy® f 1 y(Q n U»(y) \ 



In the second step we removed the explicit summation over the base points. 
A loop of length 2k without complete iteration of its contour allows for 2k 
different choices of a base point thus canceling the factor l/2k in the previous 
expression. A loop which iterates its whole contour 1(1) > 1 times allows 
only for 2k/ 1(1) different base points and a factor 1/1(1) remains (compare 
[f[7| , |l~8|1 for a more detailed discussion). We defined £ to be the set of all 
closed, connected, non back-tracking loops of arbitrary length. Note that 
each loop comes with both of its orientations and the factors from the link 
variables for the two orientations are complex conjugate to each other. 



5 



3 Expansion of the Boltzmann factor and identi- 
fication of the loop representation 



The final expression (||) for the hopping expansion performed in the last 
section is the starting point for the identification of the loop representation 
for the determinant in this section. As already outlined in the discussion 
above we will proceed by expanding the Boltzmann factor. This expansion 
converges and represents the determinant when the sum in the exponent 
converges, i.e. for m + 6{x) > as discussed above. The resulting series 
converges even absolutely and the reordering of terms which we perform 
below is justified (see e.g. ]i~9[). 

We approach the expansion by first discussing several examples which 
will be used to extract the general algebraic rules for the different terms 
along the way. Keeping track of the terms in the Taylor expansion of the 
Boltzmann factor is most easily done by using the actual contours of the 
loops associated with the different terms. So let's first look at the different 
contributions to the sum in the exponent of the Boltzmann factor. In Fig. [I] 
we show some of the terms as they appear in the exponent. 

The first two terms are two loops around a single plaquette and they 
differ only in their orientation. Note that each of these two loops (similarly 
for all loops shown) is only a representative for all loops around a single 
plaquette, the other loops differing by a translation, i.e. winding around any 
other single plaquette of the lattice. These first two loops both have four 
corners thus giving rise to a geometrical factor of (l/\/2) 4 - Certainly there 
are also the factors from taking the product of the link variables around the 
plaquette. These are however not relevant for the following discussion and 
thus are not displayed explicitly. We show only the corner-, self-intersection- 
and iteration-factors. 

Next we plot another rectangular loop, but now enclosing six plaquettes. 
Again the corner- factor is (l/\/2) 4 . The fourth term is a loop enclosing two 
plaquettes but now having a pinch structure. This loop has eight corners 
and no self-intersection such that its geometrical weight is (l/\/2) 8 - The 
next term we show is again a loop enclosing the same two plaquettes, the 
lower plaquette however is now run through with opposite orientation. This 
loop has one self-intersection, but only 6 corners, thus producing a factor of 
— (l/\/2) 6 - Note that also the weight from the link variables (not displayed 
in Fig. |l|) differs from the last contribution since the lower plaquette now 
has negative orientation and the link variables for this plaquette are complex 
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Figure 1: Graphical representation for some of the terms in the exponent 
of the Boltzmann factor. We display only the corner-, self-intersection- and 
iteration-factors. 



conjugate. Next we show a loop which runs around a single plaquette twice. 
It has eight corners and one self intersection. In addition we have to take 
into account a factor of 1/2 since the loop iterates its contour twice. The 
total factor then is — l/2(l/\/2) 8 . Finally we show another loop with a self- 
intersection. It consists of 2 loops around single plaquettes connected by a 
double line along two links. It has one self- intersection and 8 corners such 
that the geometrical factor is — (l/\/2) 8 - 

Note that in Fig. |l| we display only a few terms which are useful for the 
discussion of some examples below. It is however straightforward to depict 
arbitrary terms from the series in the exponent of (||) and to determine the 
corresponding geometrical factors. 

Let's now perform the expansion of the Boltzmann factor in (|J). We 
will proceed as follows: In the figures below we show on the left hand side 
some configuration of links as it appears in the final representation, i.e. af- 
ter the Boltzmann factor has been expanded. These contributions in their 
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final form will however be built from different loops and their products, 
emerging when expanding the Boltzmann factor. These different terms are 
shown in the middle of the figures below, together with their geometrical 
weights. The total weight associated to the graph depicted on the left hand 
side is the sum of the weights for all contributing terms. We give this num- 
ber on the right hand side, separated by an arrow. We remark that the 
expansion of the Boltzmann factor will produce contributions consisting of 
several disconnected pieces located on different positions of our lattice. Up 
to the combinatorial factors which will be discussed below, disconnected 
pieces in different areas of the lattice do not communicate with each other. 
Our examples thus concentrate on the different terms locally building up 
a contribution. The inclusion of disconnected pieces located elsewhere is a 
trivial generalization. 




Figure 2: A loop around a single plaquette. 

In Fig. [2] we show a loop which consists only of the four links around 
a single plaquette. As is obvious from Eq. @ and its geometrical repre- 
sentation in Fig. ||, there is only one contribution in the Taylor expansion 
of the Boltzmann factor giving rise to this contribution. It is the first loop 
depicted in Fig. Furthermore only the linear term in the power series 
for the exponential function produces this term since for higher powers, this 
loop would be multiplied with other loops. Thus we find only one term on 
the right hand side of Fig. ^. The total weight for this term is the same as 
we had in the exponent, namely (l/\/2) 4 - 
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Figure 3: A loop which iterates its contour twice. 

Next we consider a contribution where a loop iterates it's contour twice 
(Fig. ^). In particular we consider the simplest such term, i.e. the config- 
uration where all links around a simple plaquette are occupied twice. This 
configuration firstly can be resolved with a single connected loop. This loop 
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has eight corners and one self-intersection which produces a minus sign. In 
addition we have to take into account a factor of 1/2 since this loop com- 
pletely iterates its contour twice. The overall factor thus is — l/2(l/\/2) 8 - 
Secondly the configuration can be resolved as the square of the loop which 
runs around our plaquette only once. Each of the single loops has four cor- 
ners giving a total corner factor of (l/\/2) 8 - In addition one has a factor 1/2 
which is the factor for the quadratic term in the power series of the expo- 
nential function. Putting things together we find that the two contributions 
cancel each other. 

This type of cancellation which we have seen in the simple example of 
Fig. H is the manifestation of a general rule: Whenever a term contains a 
loop which iterates its complete contour twice the contribution gets canceled. 
The general argument goes as follows: Let l\ be the loop which iterates its 
contour twice. Assume, that it appears in a contribution which also contains 
Ti — 1 other loops li,i = 2,...n, which are located somewhere else on the 
lattice. This contribution then appears in the n-th power term of the Taylor 
series for the exponential function. The corresponding factor is: 



The 1/n! factor on the left hand side comes from the power series for the 
exponential function. Then follows a minus sign from the self-intersection 
of l\ and a factor of 1/2 since l\ iterates its contour twice. Next we have 
the corner factor for l\ followed by the factors f(h) for the other loops (cor- 
ner, self-intersection- and iteration- factors). Finally we find a combinatorial 
factor of n\ which gives the number of different ways the n terms can be 
arranged. Since h iterates its contour twice, the same contribution can also 
be obtained by squaring a loop l[ which has the same contour as 1%, but 
runs through it only once. This contribution then appears in the n + 1-th 
power term of the Taylor series for the exponential function and comes with 
a factor of 



Here the factor of —1/2 is gone, since we have no self- intersection and also 
the contour is not run through twice by a single loop. The corner factor 
stays the same since c{l\) = 2c(/^). Finally the combinatorial factor changes 
to (n + l)!/2, where the factorial gives the number of ways the n + 1 terms 
can be arranged. However, l[ appears twice, such that we have to divide 
by 2. Thus we find that (Q) is exactly the negative of (||). This proves 
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that whenever a loop iterates its contour twice it can be matched with a 
contribution which has the opposite sign and gets canceled. 

In principle one can extend this proof and show via induction that all 
contributions, where a loop completely iterates its contour more than once, 
get canceled. There is however a simpler argument. It is based on the fact, 
that not more than two lines can flow through a site. To understand this 
we have to remember that each occupied link carries with it a factor of 
1/M(x) as can be seen from the last product in (^). Thus when n loops run 
through a site xq they produce a factor of l/M(xo) n . Combined with the 
overall factor of M(xq) 2 we find M(xo) 2 ~ n . When n exceeds 2, i.e. more 
than two loops run through a site, this term is no longer polynomial. On 
the other hand one finds by direct evaluation that the expression for the 
determinant on the finite lattice is a polynomial in M(xq). Using the fact 
that for m + 8(x) > the series represents the determinant (compare the 
discussion above) we conclude, that all vertices where more than 2 lines run 
through a site have to vanish. This implies that loops which iterate their 
contour more than twice cannot contribute and for the case of the doubly 
iterated loop we have given the explicit argument for its cancellation above, 
thus showing that eventually none of the iterating loops in the exponent 
contributes. 

Before we go back to the discussion of further examples, let's bring in 
another harvest from the discussion in the last few paragraphs: The factors 
1/re! from the Taylor series of the exponential function always get canceled 
and the combinatorial factor simply is 1. This can be seen as follows: For a 
contribution which is the product of n loops one finds that only terms where 
n different loops are multiplied can contribute (if a loop could appear more 
than once we would have loops which are multiply occupied, a case which 
we have just excluded). These products come with a combinatorial factor 
of n!, which gives the number of ways the terms can be arranged. It cancels 
the factor of 1/n! from the Taylor series for the exponential function. Thus 
the algebraic factor is always equal to 1. 




Figure 4: Two nested loops with opposite orientation. 
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After this combinatorial excursion let's go back to discussing more ex- 
amples. Next we analyze a graph with two nested loops (Fig. |j). There is 
only one contribution in the expansion of the Boltzmann factor giving rise 
to this graph. It is a product of two loops, i.e. it comes from the quadratic 
term in the Taylor series for the exponential function. Thus it has a factor 
of 1/2 which is the coefficient of the quadratic term in the Taylor series. 
However, for our graph two different loops have to be multiplied, making 
our contribution a cross term which comes with a factor of 2 canceling the 
factor 1/2 from the power series. Each of the two loops has 4 corners giv- 
ing rise to a geometrical factor of (l/\/2) 8 - This is already the final factor 
associated with the contribution depicted on the left hand side of Fig. ||. 
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Figure 5: Two nested loops with equal orientation. 

In Fig. H we look at a graph which contains the same links, but the ori- 
entation of the inner loop has been reversed. In addition to the contribution 
from a product of two loops, similar to the last example, the new config- 
uration also allows for a single connected loop running through the chosen 
set of links. This term is produced by the linear term of the power series 
for the exponential function and thus simply has an algebraic factor of 1. 
The geometric factor is — (l/\/2) 8 , since the loop has eight corners and one 
self-intersection. The product contribution has the same factor as in the last 
example. The two contributions cancel each other. 

Fig. H demonstrates another general rule of our loop calculus: The final 
expression will not contain any terms where a link is occupied twice with 
the same orientation. Whenever we compute the weight for such a graph we 
find that the different terms in the expansion of the Boltzmann factor cancel 
each other. In Fig. [5| this is demonstrated in an example, but let's give the 
general argument: Whenever we have a configuration where along at least 
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one link two lines run parallel with the same orientation, then there are two 
ways to resolve this configuration with terms coming from the expansion of 
the Boltzmann factor. One contribution where the two lines do not cross and 
one with a crossing. Both contributions have the same number of corners, 
the latter one however obtains an extra minus sign due to the additional self 
intersection. The two terms thus cancel each other. We remark that at some 
other parts of our loop there also might be structures which can be resolved 
in different ways. However, already one link which is run through twice in 
the same direction leads to a vanishing of this contribution (compare also 
the discussion of completely iterated loops given above). This exclusion of 
contributions with links which are run through twice in the same direction 
can be interpreted as a manifestation of the Pauli principle which forbids 
two excitations of the same state. 




Figure 6: A configuration with a self intersection. 



Next we concentrate on how a configuration containing an intersection 
can be resolved in different ways. The left hand side of Fig. ^ shows a 
configuration which contains two plaquettes with opposite orientation which 
share a common corner. There is one possibility to run around the two 
plaquettes in a single loop. This contribution comes from the linear term 
in the Taylor series for the exponential function and has eight corners thus 
acquiring a factor of (l/\/2) 8 - Our two plaquettes can also be obtained as 
the product of two loops, i.e. a contribution appearing in the quadratic term 
of the power series. Again we have a total of eight corners giving a factor of 
(l/\/2) 8 and when adding this to the coefficient of the linear term we find 
an overall factor of (I/1/2) 6 . 

In the configuration shown in Fig. |?] we now have reverted the orien- 
tation of the lower plaquette. Certainly this configuration will still have a 
contribution from the quadratic term of the exponential function where it 
is written as a product of two loops. The geometrical factor is the same as 
in the last example. However, the linear term, where our contour appears 
as a single connected loop, has to be changed. It now has only six corners 
but one self-intersection and thus a factor of — (l/\/2) 6 . Adding this to the 
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Figure 7: Another configuration with a self intersection. 



coefficient of the product term we find an overall factor of — (l/\/2) 8 - 
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Figure 8: Loops around two plaquettes joined by a double line. 

Finally we discuss (Fig. ||) the contribution which corresponds to the 
last loop shown in Fig. |l|. It cannot be decomposed into the product of 
two loops and thus is entirely built from the linear term in the expansion 
of the exponential function. The factor — (l/\/2) 8 is simply taken over from 
the exponent. Although not very interesting from a combinatorial point of 
view, this example demonstrates another important point. In order to avoid 
overcounting of self-intersections, we have to introduce a normal ordering 
for double lines with two flows in opposite direction (the case of lines with 
equal direction has already been ruled out above). We always draw a line 
that flows from left to right above a line which flows from right to left. For 
vertical double lines we adopt the convention that downward flowing lines 
are drawn on the right hand side of upward flowing lines. 

After having analyzed several of these examples one finds that the contri- 
butions on the left hand side can be decomposed into finitely many elements 
and all contributions which appear in the expansion of the Boltzmann fac- 
tor can be built from these elements. We view these elements as tiles with 
some configurations of oriented lines drawn on these tiles (compare Fig. ^). 
The tiles can now be arranged on our rectangular lattice such that ingoing 
and outcoming flow lines on neighboring tiles match each other. To each 
tile we assign a weight Wi and the total weight of a configuration is given 
by the product of all weights for the tiles which were used in building the 
configuration. When summing over all possible ways of different tilings of 
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our lattice (placing the centers of the tiles on the sites of the lattice) we re- 
produce all possible fermion loops contributing to the fermion determinant. 
The concept of writing partition functions as a sum over tiles (vertices) was 



first introduced in [20, 21] and these models are known as vertex models 
(see also p2]|). 

The remaining problem is to assign the correct weights to our tiles. From 
the example shown in Fig. ^ we find that a simple corner has to be assigned a 
factor of l/\/2 which gives the weights W5, ...,^12 in Fig. [| Fig. |] shows that 
tiles with a straight line (wi, W4 in Fig. |^) obtain a factor of 1. The same 
factor has to be assigned when two straight lines run anti-parallel (11)44, 1^45)- 
The example discussed in Fig. || furthermore shows that antiparallel flow 
lines have to be normal ordered and our convention is to draw flow lines 
from left to right above flows from right to left and downward flows are draw 
right of upward lines. From Fig. E] we learn, that tiles with a straight line 
and a corner (1^19, W34 in Fig. ^) have a factor of l/\/2 and an additional 
minus sign if the normal ordering of the double line enforces a crossing. 
Similarly one finds that tiles with two corners (1035, 1042, 1045, 1049 in 
Fig. |) have a factor of (l/\/2) 2 = 1/2, again with an additional minus sign 
if a crossing occurs. Figs. ||,|| and the general argument given there establish 
that all tiles which contain two lines running parallel in the same direction 
must have vanishing weights (these tiles are not displayed in Fig. ^). From 
Figs. H,f?] we find that for tiles showing two crossing lines we have to assign 
a factor of —1/2 when the crossing can be resolved with a self-intersecting 
loop («?i3, wxq) and a factor of 1 when no self- intersecting resolution is 
possible (wn,wis). Finally we also have to include the empty tile (1/^43) 
which has a factor of 1. In Fig. || we give a complete list of all 49 tiles with 
non-vanishing weights. 

We remark that tiles with an odd number of incoming or outgoing lines 
of flow cannot appear, since all contributions come from products of closed 
loops and thus the net flow into a site has to vanish. Furthermore not 
more than two lines can run through a vertex as follows from the fact that 
the determinant is a polynomial in the local mass terms (see the discussion 
below Eq. (||)). The final expression for the fermion determinant reads 

det^ = n^) 2 EriK-r (t) n W^- o» 

xeA teT j=i M& yyj 

Here we introduced the set T of all admissible tilings of our lattice A with 
the tiles depicted in Fig. [9, i.e. the set of all arrangements of tiles such that 
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Figure 9: The 49 tiles (vertices) with non-vanishing weights Wi. 
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the flow lines on neighboring tiles match. By rij(t) we denote the abundance 
of tile number j in a particular tiling t £ T. The last product runs over 
all links occupied by the tiling t. For links run through in negative direc- 
tion, the link variables have to be complex conjugated. Our result is valid 
for two sets of boundary conditions: 1.) The fermions have open bound- 
ary conditions for both directions. In this case also the vertex model has 
open boundary conditions for both directions. 2.) The fermions have one 
direction with open boundary conditions and the second direction is com- 
pactified with anti-periodic boundary conditions for the fermions (cylinder 
geometry). The vertex model then has open boundary conditions for one 
direction, whereas the boundary conditions for the compact direction are 
periodic. 



Let's now briefly discuss the applicability of our method to higher di- 
mensional theories and non-abelian gauge groups. The first step will again 
be a hopping expansion leading to formula (||). As in two dimensions this 
series will converge and represent the determinant for \\H\\ < 1. Along the 
lines of the proof in [16| one can again establish this bound on the norm 
of the hopping matrix for a certain range of values for the mass m and the 
auxiliary field 9. Also the trace over the gamma matrices is known in four 
dimensions [17|. For the parameters where \ \H\\ < 1 the Boltzmann factor 
can be expanded in the same way as in 2-d, and the resulting series repre- 
sents the determinant [15|. Again we can apply the power counting trick for 
the local mass terms and show, that the number of fermion lines admissible 
at a single site is bounded. In four dimensions this bound is given by 4iV c 
where N c denotes the number of colors. For non-abelian gauge fields in ad- 
dition one has to take the trace over the product of the link variables. The 
resulting formula is the loop representation of the fermion determinant. 

Such loop representations can be seen to lead to new ways for a numerical 
evaluation of the fermion determinant by directly updating the loop configu- 
rations 1 23]. It is now possible to overcome the limitations of pseudo- fermion 
methods and simulate also odd numbers of flavors. For e.g the one-flavor lat- 
tice Schwinger model with Wilson fermions no numerical analysis exists and 
we are currently implementing algorithms based on the loop representation 
for this model |23t| . 

In this context it is interesting to note, that the loop representation for 
the fermion determinant also solves the fermion sign problem. This prob- 
lem plagues a direct evaluation of the Grassmann path integral |0] , as well 
as the Hamiltonian treatment of fermionic systems |12||. When numeri- 
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cally updating fermions, the fermion sign problem enforces an exponentially 
increasing number of Monte-Carlo configurations as the volume of the sys- 
tem is increased or its temperature is lowered |12(| . In our representation 
these permutation signs are gone and the remaining signs are due to the 
traces over the gamma matrices. These latter signs are less severe, since the 
weights for configurations with negative sign are exponentially suppressed 
due to additional factors of l/\/2 which are necessary to build loops with 
self intersections. 

Finally we argue that the loop representation could also lead to new 
methods for the simulation of lattice field theories with non- vanishing chem- 
ical potential. Although the Schwinger model with chemical potential is not 
particularly interesting per se (see e.g. [24] for a study in the continuum) 
it nevertheless poses the same technical problems as QCD with chemical 
potential. When introducing the chemical potential /j, along the lines of 
[p5| each link forward in time obtains a factor of e~ M and a factor of 
for links backward in time. Thus the symmetry between the two possible 
orientations of loops winding around the compact direction is broken and 
the determinant accquires a complex phase. This complex phase prohibits 
the application of the pseudo-fermion method which requires a non-negative 
kernel for the action of the pseudo- fermions. The loop representation allows 
for an update of single loops and the chemical potential can be included in 
a natural way independent of the number of dimensions. Also this idea is 
currently being tested using Thirring- and Gross-Neveu-type models in 
loop representation. 



4 Two limiting cases of our representation 

Our final representation @ of the fermion determinant has a relatively sim- 
ple form. It is however instructive to reduce it to limiting cases and compare 
it to known results for these instances which were obtained with different 
methods. In particular we will discuss the case of the determinant in a scalar 
background field and the case of the strongly coupled Schwinger model. 

We start with the case of the reduction to the case of the fermion deter- 
minant in a scalar background field |l^] . There the result for the fermion 
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determinant reads 



detK 



n m ^ 




(10) 



The sum now runs over the set C sa of self-avoiding, closed, non back-tracking 
loops. Self-avoiding means that the loops may neither self-intersect nor 
touch each other, they can, however, consist of several disconnected pieces. 
The loops in C sa are not oriented, i.e. each loop is counted only with one of 
its two possible orientations. c(l) again denotes the number of corners of I 
and the product in (|l(]) runs over all sites y visited by the loops. 

Finding a general proof which directly shows the equivalence of the right 
hand side of (|9|) at U^x) = 1 with the right hand side of ([R]) turns out 
to be a rather hard combinatorial problem. Instead we here discuss an 
example which illustrates how the two representations are related. The 
main difference is that the representation (||) for the determinant where 
also a gauge field is coupled has to work with oriented loops, since the 
link variables allow to distinguish between forward and backward hopping 
(complex conjugation in the kernel (|l|)). Loops in the determinant for a 
scalar background field give the same contribution when their orientation is 
reverted, and thus it is possible to find the representation (|n]) in terms of 
non-oriented loops. 
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+ 



+ 



Figure 10: The contributions to a given configuration as they appear in 
the scalar formula (u( 



Let's now look at a simple configuration which we depict in Fig. |10J. In 
the figure we show how this term is generated in formula (|i~0|), i.e. how it 
can be written as the product of two self-avoiding loops. There are three 
ways of doing so, each of them is a cross term in the square, such that we 
collect an overall factor of 2. The first term gives rise to 8 corners, while 
the other two configurations each have a total of 12 corners. We thus find a 
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Figure 11: The same contribution as in Fig. now being built from the 
oriented loops of formula (||). 
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Next we generate the same term, now using the tiling formula (||). Again 
we find 3 topologically different configurations which we show in Fig. [Tl]. 
In addition these three configurations come in pairs related to the plotted 
terms by reverting all arrows. This gives again an overall factor of 2. By 
collecting the factors for the vertices as given in Fig. we find 



/ 1 \8 / 1 \8 1 1 
I \ y y + 1 



, _ X 1 X 1 + I —= I X - X , _ , 

v^2^ 2 2 v^y 



+ ( —) x - X 



(- 



1 1 

- x - 

2 2 



3 
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Thus both formulas (jpj) and ( |10|) give the same factor. Similarly one can 
show the equivalence of (^) and (^) for more complicated examples. 



Let's now proceed to the case of the strongly coupled Schwinger model. 
All scalar field variables are set to zero, 9{x) = 0, and we write the U(l) 
link variables now explicitly as 



Ap(x) G (-vr,+vr] 



The partition function of the lattice Schwinger model is then given by 



Z n 



jj dAjAx)_ e -^s gauge [A] dH K 



2tt 



(11) 



At strong (= infinite) coupling g, the Boltzmann factor for the gauge field 
equals 1. Only the determinant appears in the integrand. Each t in the sum 
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over all tilings T gives rise to links being empty, singly occupied links and 
doubly occupied links with link variables complex conjugate to each other. 
These cases correspond to the integrals 



- dA^x) = i r dA^{x) c±iAu[x) 

-7T 1~K J-tt 2lT 



0, 



71 dA^(x) +iAu (x)-iA u (x) 



2tt 



1. 



The second integral implies that whenever a configuration of tiles contains 
a link which is occupied only once, the whole configuration vanishes when 
integrating the gauge field. Thus only configurations with either empty or 
doubly occupied links (with opposite orientation) contribute in the sum over 
all tilings, i.e. only the tiles 43, ...,49 in Fig. ^ contribute. The strong cou- 
pling limit thus reduces the lattice Schwinger model to a seven vertex model 



with vertices and weights as depicted in Fig. 12, where we have included 
the overall factor Yl x (2 + m) 2 into the monomer weight w\. Note that after 
integrating the gauge fields in the strong coupling limit the orientation of 
the link variables gets lost. Our loops do no longer have an orientation and 



it is sufficient to have single lines on the tiles in Fig. 12 



Wj= (2+m) w 2 =l 



w 3 =l 



r 



n 



w=l/2 w=l/2 

4 5 



w=l/2 w=l/2 

6 7 



Figure 12: The 7 vertices occuring in the strong coupling limit and their 
weights. 

The partition function in the strong coupling limit reads 



Jg=oo 



E 

ter 7 j=i 



n( 



(12) 



Here Tj denotes the set of admissible tilings obtained from the seven tiles 
depicted in Fig. 12. Formula (12) is equivalent to the partition function for 
self- avoiding loops with corner activity 1/2 and monomer weight (2 + m) 2 . 



The representation (Il2 
found by Salmhofer 



for the strongly coupled Schwinger model was first 
. Since || uses different techniques (expansion of 
the Boltzmann factor and local saturation of the Grassmann integral) the 
derivation of this result as a limiting case of our general formula provides a 
second non-trivial test of Eq. (pi). 
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It is interesting to note that our representation (^) allows to generalize 
the vertex model representation for the strongly coupled case to a model with 
N flavors. For N flavors the integration of the fermions gives the iV-th power 
of the fermion determinant. This simply corresponds to superimposing N 
different 49-vertex models with loops in N different colors onto each other. 
Integrating the gauge fields in the strong coupling limit discards all vertices 
where at least one link is only singly occupied. For the case of the strongly 
coupled two-flavor Schwinger model this leads to a 163- vertex model p6| . 

We remark that the representation of the strongly coupled lattice Schwin- 
ger model as a seven vertex model demonstrates the power of loop represen- 
tations when performing numerical simulations. The model ( |l2|) has been 



studied numerically in p7|, 28] and high precision data on the phase dia- 



gram has been obtained. For vertex models furthermore the use of cluster 
algorithms might become possible which could further increase the quality 
of numerical data [29]. 



5 Integrating out the gauge fields at arbitrary cou- 
pling 

We now integrate the gauge fields at arbitrary coupling for the model with 
open boundary conditions. The basic idea is to fill the loops with plaquette 
variables. We remark that when choosing our lattice to be the above dis- 
cussed cylinder we find topologically non-trivial loops. For these loops the 
gauge field integral, although still solvable in closed form, is more involved 
and for convenience we here discuss only the case of a rectangular lattice 
with open boundary conditions. 

We now assign plaquette variables V(p) to each plaquette p, which are 
given by the product of the 4 link variables associated to the plaquette p 
when running through it in mathematically positive direction, 

v(p) = n w ■ ( i3 ) 

(x,n)ep 

To each plaquette p we furthermore assign an occupation number J p which 
gives the power of V(p) for the plaquette p. Negative J p is used when we 
assign the complex conjugate plaquette variable V(p)* with exponent |J p |. 
When two neighboring plaquettes have the same occupation numbers, the 
link variable shared by the two plaquettes cancels. When the occupation 
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numbers differ, the shared link variable remains. Note that for neighboring 
plaquettes the occupation numbers can only differ by dbl, since we can never 



have more than one unit of flux on a link. Fig. 12 gives an example of a 



configuration of loops and how its occupation numbers are assigned. 
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Figure 13: A contribution to the determinant and the occupation numbers 
for the plaquettes inside the loops. 

For a given configuration of occupation numbers J p the gauge field inte- 
gral reads, 



Wp) = J n du > 



(x) w n ( n Mv) 



(14) 



In our two-dimensional model this integral can easily be solved by choosing 
a simple gauge, i.e. a maximal tree of link variables set to 1, which factorizes 
the gauge field integral. The gauge field integral is simply a product over 
all plaquettes with the factors determined by the occupation numbers, 



(15) 



where the function F g (n) , n S 7L is given by (we quote the results for Wilson 
action as well as for Villain action), 



F g {n) 



e ^ l\ n \ (jz) (Wilson action) , 



(16) 



9_ e -V»' 



2t 



(Villain action) . 
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Here I n denotes the modified Bessel function. For the partition function of 
the lattice Schwinger model we thus find, 

r i 49 

Z 9 = / II dU^e-^^detK = £ l^w^ IIWpW) ■ 

(x,m) *er j=i p 

(17) 

The last term denotes the product over all plaquette factors depending on 
the occupation numbers J p (t) for a given tiling t. We remark that we now 
have included the factors of (2 + to) into the new weights vjj. They are 
obtained by multiplying the old weight W43 by (2 + to) 2 , by multiplying 
weights wi, ...,wi2 with 2+to and setting Wj = Wj for the remaining weights. 

Our expression ( |i~7| ) represents the partition function as a sum over dis- 
crete variables. This expression could lead to the possibility of applying 
cluster algorithms [29| which might considerably improve the quality of nu- 
merical data. 



6 Observables 

In this section we discuss the loop representations for pure gauge field ob- 
servables such as the Wilson loop and for local fermionic observables in 
particular scalars and vectors. 

Let's start with the simpler observables which are the ones only con- 
taining gauge fields. The prototype of such an observable is the Wilson 
loop 

w c (u) = n u,{x) . 

Here C is a closed contour and without restriction of generality we assume 
that it is oriented in mathematically positive direction. As for the fermion 
loops, we can create the link variables on C by filling the whole contour with 
plaquette variables as defined in ( |l3|) . This requires an occupation number of 
+1 for all plaquettes inside C. This occupation number can simply be added 
to the ocupation numbers J p from the fermion loops. When integrating out 
the gauge fields we find for the vacuum expectation value of the Wilson loop 
in the lattice Schwinger model (compare (|l7|)) 

(Wc(U)) g = ttE I[(^) nj(t) I[F g (j p (t) + Xp (C)) , (18) 
^9 teT j=i p 



23 



where we defined Xp(Q to be the characteristic function of the area inside 
the chosen contour C with x«(C) = 1 for plaquettes p in the interior of C and 
Xp(C) = for plaquettes outside. It is obvious that this formula can easily 
be generalized to more involved gauge invariant observables for the gauge 
field. 



Only slightly more complicated is the loop representation for local fermi- 
onic observables, in particular scalars and vectors. Here we can make use of 
the fact, that our formalism gives the fermion determinant, i.e. the fermionic 
path integral, in an arbitrary gauge field and in an arbitrary local scalar 
field 9{x). Thus by differentiating the determinant with respect to these 
external fields we can generate n-point functions for the fermionic bilinears 
in the action, i.e. scalars and vectors. By writing the link variables as 
Ufj,{x) = exp(iA^(x)) with A^ a (x) = —A^(x — fi), the kernel K for the 
fermion action reads (compare (Q)) 



K(x,y) = 2 + m + 8(x) 



±2 



u=±l 



Thus we obtain for the expectation value (just the fermion integration) of 
scalars and vectors 
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detK 



The current j u (z) is the Noether current corresponding to the gauge sym- 
metry of our action 



il)(z)<T v U v (z)rl){z + v)+ 4>(z + v)<j v U v (z)*i){z) 



il){z)U„{z)^{z + v) - i>(z + u)U u {z)*^{z) 



(19) 
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In the continuum limit the first term turns into the vector-current i/j(z)'y v il}(z) 
while the second term is of order 0(a) and vanishes. 

Now we act with the differential operators d/d9{v) and d/dA u (z) on our 
loop expression (|9|) . When making explicit all external fields the formula (^) 
reads (now again using the original vertex weights Wi from Fig. |9| and not 
the rescaled vertices Wi used in ([FT 



49 P ±iA„{y) 

In the numerator of the last term which collects the link factors, the plus 
sign in the exponent has to be taken for links (y, fi) which are run through 
in positive direction for the configuration t and the minus sign is taken for 
links occupied in negative direction. When acting with d/d9(v) for some 
fixed site v on the expression (^) we have to distinguish three different cases. 
When two fermion lines run through v, all terms containing 9(v) cancel each 
other and the differentiation simply gives zero. When one fermion line runs 
through v this contribution is linear in 2 + m + 9{v) and after differentiation 
there remains a factor of 1. Finally when the site v is empty the factor 
in (p0|) is [2 + m + 9{v)] 2 and the action of the differential operator then 
gives a weight of 2 [2 + m + 9(v)]. The results for these three cases can be 
denoted conveniently by inserting the following factor in the sum over the 
configurations t S T (we have now set the scalar field 9 to zero everywhere 
which gives us the expression for the Schwinger model) 

where we defined si n \t) to be equal to 1 when the site v is visited exactly 
n-times by the fermion loops of the configuration t and zero otherwise. The 
first factor in (|2l] ) produces a zero whenever 2 fermion lines run through v. 
The other two terms correct the weight at v for the cases of a single fermion 
line and for v being empty. 

The result for the differential operator id/dA u (z) is even simpler to de- 
note. From (|20| ) it is clear that for configurations t where the link (z, v) 
is empty, or configurations where this link is doubly occupied and the link 
variables cancel each other due to their opposite orientation, this contri- 
bution does not depend on A v {z) and is annihilated by id / dA u (z). For 
configurations where the link (z, v) is run through with positive orientation, 
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we obtain after applying id/dA u {z) a factor of — 1 and a factor of +1 for 
negative orientation. Thus the differential operator simply measures the 
negative of the net flow of fermion lines through the link (z, v). By defining 



+ 1 when t occupies (z,u) in positive direction 
L( z ^ v ){t) = < — 1 when t occupies (z,u) in negative direction 
when t has net flux zero through (z, v) 



(22) 



we can take into account the vectors by simply including —Lr z>v \(t) in the 
sum. Inserting this factor, together with (21) into the sum over all configu- 
rations we obtain for the fermionic n-point functions in the lattice Schwinger 
model (the gauge fields were integrated out as above, 9{x) is set to 0, and 
we use again the rescaled weights w-i as defined below equation fll7|)) 
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\{L {zuUl) (t) .(23) 
1=1 



After getting used to it, this expression turns out to be very powerful and 
also very simple from a numerical point of view. For computing n-point 
functions of scalars one merely has to locally change the weights on the 
sites supporting the scalars, and for vectors only the net flux through the 
links associated with the currents has to be counted. This is certainly much 
simpler than inverting the fermion matrix as is necessary in the standard 
numerical approach. 

It is instructive to play a little bit with formula (|23|), so let's e.g. show 
that {ju(z)) g = 0. As discussed above only configurations where the link 
(z, v) is singly occupied can contribute. These configurations, however, come 
in pairs related by reversing the orientation of all fermion lines. The two 
contributions have the same weight from the gauge fields (J p — > —J p ) but 
are counted with opposite sign in (p3|) and thus the one-point function of 
the current vanishes as it should. 

Furthermore it is interesting to analyze some variables in the strong cou- 
pling limit, in particular the chiral condensate (ipip(v)) g=00 . As discussed 
above, in the strong coupling limit all configurations with singly occupied 
links are suppressed. Thus, according to (|2l"|), the chiral condensate can 
obtain contributions only from the empty vertex. Using translational in- 
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variance of {ij)ip{v)) g=00 on the infinite lattice, one finds (compare Eq. 12) 

2 



g=oo 



2 + m A- 



lim 



1 



EIIk v '" 



ni(t) 
|A| ■ 



(24) 



Thus in the strong coupling limit the chiral condensate is proportional to 
the density of vertex 1 of Fig. 12, i.e. the empty tile. Similarly it is possible 
to obtain simple expressions for other fermionic observables in the strong 
coupling limit. 



7 Discussion 

In this article we have derived the loop representation for the fermion deter- 
minant of 2-D Wilson lattice fermions coupled to U(l) gauge fields, i.e. the 
fermion determinant of the lattice Schwinger model. The construction is 
based on first using standard hopping expansion which already expresses 
the determinant in terms of loops which are the natural (gauge invariant) 
variables for the fermion determinant. Standard hopping expansion, how- 
ever, only gives a loop expansion for the free energy, and thus to obtain the 
desired result, the Boltzmann factor was expanded. An essential ingredient 
for this expansion was the introduction of a locally varying mass term. For 
our problem it was shown that the determinant is a polynomial in the local 
mass terms with a degree of two. This considerably restricts the number 
of terms which can appear when expanding the Boltzmann factor. For the 
remaining terms it was shown that they can be obtained by summing over 
local pieces of loops, i.e. over the vertices given in Fig. ^. 

The gauge fields were integrated out and the result was expressed in 
terms of occupation numbers for the plaquettes enclosed in the fermion 
loops. The whole formalism generalizes to the case of non-vanishing chemi- 
cal potential in a natural way. Gauge field variables can be treated on the 
same footing as the fermion loops, i.e. by assigning occupation numbers for 
the plaquettes followed by integrating the link variables. Local fermionic 
observables were generated by differentiating the determinant with respect 
to the external fields. In the loop picture their vacuum expectation values 
simply correspond to counting the net flow of fermion loops through the 
sites (for scalars) and links (vectors) of the lattice. 

Certainly the lattice Schwinger model is a relatively simple system, but 
the essential ingredients of our construction can be carried over to more 
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realistic models. In particular, a hopping expansion of the fermion determi- 
nant for Wilson fermions is straightforward also for higher dimensions and 
non-abelian gauge groups. Furthermore the necessary traces over the 7- 
matrices can be computed in closed form also in 4 dimensions jL7|] . It is also 
possible to again apply the trick of introducing a locally varying mass term 
and, in the region of convergence, use it to rule out infinitely many terms in 
the expansion of the Boltzmann factor. The number of allowed loops run- 
ning through links and sites will depend on the dimension of the problem 
and the number of colors. We expect that also in higher dimensions it is 
possible to decompose the loops into local elements. For non-abelian mod- 
els, these elements will, however, be more involved, since also for the gauge 
fields one now has to compute traces. The inclusion of chemical potential 
can be treated in exactly the same way as we have done for the Schwinger 
model. Also local fermionic observables will again correspond to counting 
the net flow of fermion lines through sites and links. The treatment of the 
gauge field integral, however, becomes more complicated. For abelian fields 
in higher dimensions it is still possible to build the links on our loops by 
filling surfaces ||, and a subsequent duality transformation (see e.g. f30|| ) 
produces simple Gaussian weights. For non-abelian gauge groups no such 
techniques are known and probably one would have to rely on perturbative 
calculations. 

We remark that also fermionic models such as the Thirring model or 
the Gross-Neveu model are accessible to the techniques developed here. 
The four-fermi terms are then generated by coupling an auxiliary field to 
fermionic bilinears and these auxiliary fields are subsequently integrated 
with Gaussian measures. The auxiliary fields necessary for generating the 
above mentioned models are abelian fields and also their Gaussian integrals 
are straightforward to solve (see e.g. for a study of this approach to the 
Gross-Neveu model in 2D). Thus for purely fermionic models the tools for 
finding a loop representation are readily available. It is hoped that such a 
representation opens new ways for numerical simulations of these systems, 
in particular for the case of non-vanishing chemical potential. 
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